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ABSTRACT 

Reinforced Concrete (RC) has been widely used in construction of infrastructures for many decades. The cracking 
behavior in concrete is crucial due to the harmful effects on structural performance such as serviceability and durability 
requirements. In general, in loading such structures until failure, tensile cracks develop at the initial stages of loading, 
while shear cracks dominate later. Therefore, monitoring the cracking modes is of paramount importance as it can lead 
to the prediction of the structural performance. In the past two decades, significant efforts have been made toward the 
development of automated structural health monitoring (SHM) systems. Among them, a technique that shows promises 
for monitoring RC structures is the acoustic emission (AE). This paper introduces a novel probabilistic approach based 
on Gaussian Mixture Modeling (GMM) to classify AE signals related to each crack mode. The system provides an early 
warning by recognizing nucleation of numerous critical shear cracks. The algorithm is validated through an experimental 
study on a full-scale reinforced concrete shear wall subjected to a reversed cyclic loading. A modified conventional 
classification scheme and a new criterion for crack classification are also proposed. 

Key words: Acoustic Emission, Crack mode classification, Cluster Analysis, Gaussian Mixture Modeling, Probabilistic 
damage diagnosis, Reinforced concrete. 



1. INTRODUCTION 

Structures made of reinforced concrete (RC) are ubiquitous in modern infrastructure systems. These structures are 
subjected to deterioration due to aging, increased load, and natural hazards. Conventionally, the assessment of civil 
structures relies on visual inspection (VI) , which in some cases is accompanied with some form of ground-based 
nondestructive evaluation (NDE) method. Nonetheless, severe degradation could go undetected. Notable examples of 
such unanticipated failures that have exposed shortcomings in current assessment methodologies can be found 
worldwide, such as the viaduct of the Belfast-Dublin railway line collapsed in 2009 and the Christchurch CBD buildings 
collapsed in 201 1 2 ' 3 . Structural health monitoring is a field that has received significant interest over the past few years 
and has led to the development of a wide variety of systems that can overcome the qualitative nature of Vl-based 
condition assessment 4 ~ 13 . In particular a technique that shows promises for monitoring crack propagation in RC 
structures is Acoustic Emission (AE). The past few years have seen a surge of research aim at determining functional or 
correlation interdependencies between crack characteristics and AE features 4 ' 14 , including Z>-value and Sifted Z>-value 
(SZ>) analysis 5 ' 8 ' 15 and AE energy 16 . In this paper a new approach based on a Gaussian Mixture Modeling (GMM) is 
presented. GMM is an unsupervised distribution-based classification technique 17 that has been successfully used in 
many fields, including sound recognition ' , image processing , dynamical systems and tracking , and text 
recognition 22 ; however the application of this method in AE -based structural health monitoring (SHM) systems has not 
been investigated. Such system, coupled with remote facilities 23 , could be employed to continuously monitor complex 
structures and alleviate hundreds of inspection obstacles stemming from accessibility restrictions, complex geometries, 
and the location of hidden damage. The rest of this paper is organized as follow: the next section describes the 
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conventional crack classification in RC structures and its drawbacks. The novel approach and the theory behind GMM 
are then given in section 3. Experimental setup is explained in section 4 followed by the results of structural response 
and AE monitoring in section 5. Finally conclusions are summarized in section 6. 

2. CONVENTIONAL CLASSIFICATION 

Extracting features from AE signals is usually referred to as parameter-based technique. The Japan Construction and 
Material Standard (JCMS) 24 proposed a technique to monitor the crack propagation in concrete structures based on two 
AE parameters, namely "Average Frequency (AF)" and "RA" value. AF is defined as the ratio of counts to duration and 
RA is the proportion of rise time to peak amplitude. Figure 1 illustrates this conventional crack classification approach 
using these two features. In particular, a straight line recommended line to separate tensile signals from shear signals 5 ' 7 ' 
"~ 28 . The JCMS technique is based on the observation that, in general tensile cracks generate AE waveforms with short 
rise -time and high frequency, whereas shear cracks usually result in longer rise-time with lower frequency and longer 
rise time 29 . Several experimental studies demonstrated the effectiveness of these two features 8 ' 25 ' 27 ~ 29 ; however, it 
should be mentioned that a defined criterion on the proportion of RA value and average frequency (characteristics of the 
separator line in Figure 1) has not been established yet 5 ' 7 ' 27 . In this paper, a new probabilistic approach based on a 
Gaussian Mixture Modeling (GMM) is proposed to classify AE signals. 




RA value 

Figure 1. Conventional crack classification in JCMS-IIIB5706 code 24 

3. GAUSSIAN MIXTURE MODELING 

A Gaussian Mixture Model is a parametric probability density function composed of weighted sum of M component 
Gaussian densities. For a Z)-dimensional observation or training vector X , the mixture density is defined as 18 , 
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where &>, , i=l, ..., M, are the mixture weights and 5V,- (x | // ; ,Z,- ) , i = 1, ..., M, are the component Gaussian (Normal) 
densities. Each component density is a Z)-variate Gaussian function of the form, 
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The parameters of a complete Gaussian mixture model are the mean vectors, covariance matrices and mixture weights 
from all component densities, represented together as, 
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The next step is to find the best estimate for X using Maximum Likelihood (ML) 18 ' 30 . For a sequence of T training 
vectors X = {x x , ..., x T } , the GMM likelihood can be written as, 



p(X\A) = Ylp(x t \A) 



(4) 



t=\ 



ML parameters can be obtained iteratively by Expectation Maximization (EM) algorithm 30 ' . This algorithm consists of 
two steps, Expectation and Maximization, which guarantees a monotonic increase in the model's likelihood value . 
Expectation step produce the a posteriori probability for component i, which is defined as the probability of state i given 



the observation x, and the k' re-estimated model X': 
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The maximization step returns the distribution parameter with these components 18 
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The use of a GMM is motivated by its capability to model some underlying set of hidden classes 18 ' 32 ' 33 . To implement a 
system for crack mode classification in concrete structures, the features (or measurement) is a 2-D vector (i.e., 
x = (RA, AF) ) and sequence of T training vectors is, 

X = {x x =(^,^),...,x r =(RA t ,AF t ),...,x T =(RA T ,AF T )} (9) 

and the latent or hidden classes are ;-{l,2}, that is, tensile and shear mode, respectively. 

4. EXPERIMENTAL SETUP 

The specimen was a large scale RC shear wall with a height to width ratio of 0.94, designed based on ACI 318-08 34 . 
Details about the specimen and AE setup are summarized in Table 1 . 

Table 1 . Experimental details 



RC shear wall 




Acoustic emission 


Height [mm] 


3300 


AE system 


PAC Micro-II 


Width [mm] 


3050 


AE sensor 


PAC R15a 


Thickness [mm] 


200 


Number of sensors 


8 


Reinforcement Ratio [%] 


0.67 


Preamplification [dB] 


40 


Concrete compressive strength f c [MPa] 


24.8 


Signal bandpass filter [kHz] 


20-400 


Reinforcing bars yield strength [MPa] 


464 


Trigger threshold [dB] 


35 


Reinforcing bars ultimate strength [MPa] 


708 


Sampling rate [MHz] 


1 



The main components of the AE system included an eight-channel high-speed data acquisition board (Physical Acoustics 
Corporation Micro-II PAC) and dedicated software for signal processing and storage (AEwin). AE sensors were attached 
to one face of the wall using hot glue. The experimental setup and sensor layout for loading test are shown in Figure 2. 
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Lateral load was applied to the specimen by two horizontally inclined high force capacity actuators. The specimen was 
subjected to a displacement controlled quasi-static reversed cyclic loading, see Table 2. Since very low AE activity was 
recorded prior to LS2, data in LSO and LSI was neglected in the post-processing analyses. 



Table 2. Load protocol 
























Load step 


LSO 


LSI 


LS2 


LS3 


LS4 


LS5 


LS6 


LS7 


LS8 


LS9 


LS10 


Number of Cycles 
Loading Time [sec] 
Proposed displacement [mm] 
Actual displacement [mm] 


3 

2.5 
1 
0.8 


2 
3 

1.9 
1.6 


2 
6 

3.8 
3.2 


2 
9 

5.7 
5.1 


2 
12 
7.6 
6.8 


2 
16 
10.2 
9.2 


2 

20 
12.7 
11.6 


2 

30 
19.1 
17.9 


2 

40 
25.4 
24.4 


2 

60 
38.1 
37.1 


2 

80 

50.8 

47.5 






Figure 2. Laboratory test setup and AE sensor layout 



5. RESULTS 



5.1 Structural response 



The force-displacement hysteresis loops and the corresponding backbone curves are illustrated in Figure 3. The 
specimen had a nonlinear response in LS7 and reached its ultimate strength in LS9. In the initial load steps the cracks 
started from tensile mode and then they were followed by diagonal shear cracks in the intermediate load steps. The 
diagonal cracks became dominant at the final load steps. A comprehensive report on the crack propagation pattern of this 
wall is available in 5 for interested readers. 



Proc. of SPIE Vol. 8692 86920B-4 
Downloaded From: http://proceedings.spiedigitallibrary.org/ on 04/20/2013 Terms of Use: http://spiedl.org/terms 



i, 

a 

o 

Oh 




-1500 L 



-60 -40 -20 20 40 
Displacement [mm] 

5 



60 80 



Figure 3. Force-displacement hysteresis loops and backbone curves 

5.2 GMM clustering outcome 

This section presents the results of the crack classification approach using GMM. The observation dataset X, consisting 
of x t =(RA i ,AF i ) is X = {(RA l ,AF 1 ),(RA 2 ,AF 2 ),...,(RA n ,AF n )}, Xi m 2 . As an example, X for LS4 is shown in Figure 

4. It can be observed that, using the JCMS classification scheme, establishing the best line to cluster the AE signals in 
two classes is not an easy task. 
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Figure 4. Typical pattern of Average Frequency vs. RA value and ambiguity of best separator line based on JCMS code 

To overcome this limitation, the proposed clustering analysis was applied to partition the data set X in two subsets (i.e., i 
=1, 2), named cluster 1 (i.e., tensile cracks) and cluster 2 (i.e., shear cracks). The corresponding contour of Gaussian 
mixture density function at LS4 is illustrated in Figure 5. The following observation can be made: 1) the two hidden 
classes have a mixed Gaussian density function with a larger weight for the tensile class. This is the optimum equivalent 
of data set distribution with Gaussian mixture density, guaranteed by the GMM algorithm; 2) although we will propose a 
new criterion for clustering the two hidden classes, a line could be still used by defining the segment bisector of a virtual 
line that connects the means of each class. The example of this bisector is superimposed on Figure 5 by a dashed line. 
GMM analysis was carried out for each load step separately. Figure 6 illustrates the Gaussian mixture density functions 
at load steps 2, 5, 7, and 9. At LS2 as an initial load step, tensile cracks are clearly dominant. By increasing the applied 
load at intermediate load steps, e.g., LS5, shear signals start to increase as more diagonal cracks appear 5 . At LS7 almost 
all observations are clustered in shear mode as it is shown in Figure 6c with almost no secondary summit for Gaussian 
Mixture. At the final load steps (e.g., LS9), shear cluster begin to become fatter as shown by more topographical lines on 
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Figure 6d. These results confirmed three stages of crack mode propagation: dominancy of tensile cracks at the initial 
load steps, transition stage from tensile to shear in intermediate load steps, and dominancy of shear cracks at final load 
steps. 
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RA [usec/dB] 
Figure 5. Mixed Gaussian density function of tensile and shear cracks at LS4 and the conventional separator line 
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Figure 6. Gaussian Mixture Modeling of clustering parameters for load steps 2, 3 (dominancy of tensile cracks), 5 and 6 
(Transition stage from tensile to shear), 7 and 9 (dominancy of shear cracks) 
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To improve the classification results, a new criterion is proposed to keep the data points with higher probability into 
the tensile or shear clusters and associate the rest of date set in another cluster namely "Mixed". At this aim, the 
likelihood ratio test was applied. The discriminate function is based on the log given as, 



like(X) = p(A t | X) = — i — k,FX u = exp(loglike(X)) 
p{X) 

loglike{X)=-^{X- V iJY Jk \x-^) + \n{p{l t ))-U n \Y Jt \-^\n2K 



y/ = -2 In 



like(X) for \stModel(Tensile) 
like(X)for 2ndModel(Shear) 



(10) 

(11) 

(12) 



where y/ is the likelihood ratio, fi k , V , and p(A k j are the mean, covariance, and weight of the components. The next 

step is to assign an arbitrary ratio. The value of y/j— -1.4 is selected to separate the date in which the likelihood of tensile 
is twice the likelihood of shear (-21n(2)=-1.4); likewise, the value of ^=1.4 means the likelihood of shear is twice the 
likelihood of tensile. Figure 7 illustrates how the likelihood ratio test segregates the dataset in LS5. The diamonds (0) are 
tensile signals with higher probability. The gray dots (.) are "Mixed" mode signals since they include the data with less 
probability of being either shear or tensile. The circles (o) are the shear signals with higher probability. Therefore, the 
new criterion for boundaries of clusters is defined based on a likelihood ratio value. 
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Figure 7. Data selection with higher probability in each cluster for LS5 

Figure 8 shows the number of AE signals associated to each cluster versus drift ratio. Below 0.4% drift ratio where the 
wall was completely elastic (LS5) the majority of AE signals were originated from the nucleation of tensile cracks while 
after 0.4% drift at LS6, shear signals begin to surpass tensile mode. These observation were confirmed by visual 
inspection of crack patterns in each load step 5 ' 7 . Finally, beyond 0.6% drift at load steps 7 to 9, increased number of 
shear cracks was observed. Another worthy observation can be obtained by comparing these results to the force- 
displacement hysteresis loops (Figure 3). Onset of yielding occurred in LS7; interestingly the GMM identified high AE 
activity associated to shear mode in this load step (LS7). The following increase in reported tensile cracks at LS8 can be 
due to the dominant mechanism of brittle fracture in compression that is basically Mode I cracking, the same as in 
tension. Under compressive loads, brittle materials fail by the formation of tensile microcracks at microdefects such as 
cavities and grain interfaces 35 . Table 3 summarizes the ratio of the number of data in each cluster to total data set in each 
load step. Based on these observations, a Structural Health Monitoring (SHM) system can be designed for similar 
reinforced concrete shear walls by designating a threshold on the percent of shear mode cracks. The proposed ratio for 
dominancy of shear class over the other two is 60%. Once the ratio passed that threshold, an early alarm can be triggered 
as a notice of yielding. Thereafter, imminent failure should be prevented by selecting an appropriate retrofitting scenario. 
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In addition, this classification scheme allows visualizing how much the classes are intersected. For instance, in LS2, LS7 
and LS9 the dominancy of one class of signals is comprehensible while during transition stages (LS5 and LS6) these two 
classes are highly intersected. Based on these observations we propose a new classification scheme, as shown in Figure 
9. The intersection area can be considered as mixed crack mode. 
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Figure 9. The new classification scheme 



Table 3. Percentage of AE activities for each group 



Load step 


Tensile 


Shear 


Mixed 


2 


59% 


33% 


8% 


3 


44% 


34% 


22% 


4 


56% 


34% 


10% 


5 


40% 


38% 


22% 


6 


32% 


40% 


29% 


7 


20% 


74% warning 


6% 


8 


46% 


38% 


16% 


9 


21% 


63% 


16% 



6. CONCLUSION 

This paper presented a new classification approach based on a Gaussian Mixture Modeling (GMM) of acoustic emission 
for monitoring crack mode (i.e., tensile or shear) in reinforced concrete structures. A full scale RC shear wall was tested 
to evaluate the robustness and applicability of the proposed system. Acoustic emission waveforms were recorded and 
two features, RA and Average Frequency (AF), were extracted. These data were used as input to the GMM algorithm to 
cluster the two hidden classes of shear and tensile. The results shown that three stages in the crack mode progression can 
be identified: dominancy of tensile cracks at the initial load step, transition stage from tensile to shear, and dominancy of 
shear cracks at final load steps. Conventional visual inspection confirmed the outcomes of the proposed approach. 
Further studies should investigate the effects of different aggregate mixture, compressive strength, reinforcement ratio, 
dimension, AE sensor characteristics, and loading directions (bending, tension, compression) to verify the independency 
or dependency of the proposed method on each of these factors. 
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